Fictitious time wave packet dynamics: I. Nondispersive wave packets in the quantum 

Coulomb problem 
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Nondispersive wave packets in a fictitious time variable are calculated analytically for the field- 
free hydrogen atom. As is well known by means of the Kustaanheimo-Stiefel transformation the 
Coulomb problem can be converted into that of a four-dimensional harmonic oscillator, subject to a 
constraint. This regularization makes use of a fictitious time variable, but arbitrary Gaussian wave 
packets in that time variable in general violate that constraint. The set of "restricted Gaussian 
wave packets" consistent with the constraint is constructed and shown to provide a complete basis 
for the expansion of states in the original three-dimensional coordinate space. Using that expansion 
arbitrary localized Gaussian wave packets of the hydrogen atom can be propagated analytically, and 
exhibit a nondispersive periodic behavior as functions of the fictitious time. Restricted wave packets 
with and without well defined angular momentum quantum numbers are constructed. They will be 
used as trial functions in time-dependent variational computations for the hydrogen atom in static 
external fields in the subsequent paper [T. Fabcic et al., submitted]. 

PACS numbers: 32.80.Ee, 31.15.-p, 04.30.Nk, 04.20.Jb 



I. INTRODUCTION 

Wave packets play an important role in the description 
of atoms, e.g., for the understanding of ionization pro- 
cesses in microwave experiments [l|, Q or in experiments 
with short laser pulses [1, Q. Theoretically, the wave 
packet propagation can be calculated by exact quantum 
computations Q or approximately with, e.g., semiclassi- 
cal lu\ or variational >7\ techniques. Contrary to the har- 
monic oscillator, where Gaussian wave packets and coher- 
ent states [§| can easily be described analytically, the evo- 
lution of arbitrary Rydberg wave packets is nontrivial al- 
ready in the pure Coulomb problem where spreading and 
revival phenomena are observed in the long-time com- 
puter simulation of a quantum wave packet Q. Coher- 
ent states for the hydrogen atom have been constructed 
by Klauder [lo'| and by Majumdar and Sharatchandra 
[lH, however, Bellomo and Stroud [13, [13] have shown 
that these states do not move quasiclassically but spread 
rapidly over the Keplerian orbit. Dispersion is a gen- 
eral property of Rydberg wave packets with the excep- 
tion of nondispersive electronic wave packets existing in 
periodically driven atoms such as the hydrogen atom in 
microwave fields [lH . 

The success of applying variational methods to wave 
packet propagation crucially depends on the choice of 
the trial function. Gaussian wave packets (GWPs) 
are certainly well suited for smooth and nearly har- 
monic potentials [H, [13] ■ The Coulomb potential is 
not a promising candidate for successfully propagating 
GWPs directly. Nevertheless, the GWP method based 
on the local harmonic approximation has been applied 
in one dimension to the singular Coulomb potential 
[Tsl . [TqI . [2^ . In the three-dimensional space a regulariza- 
tion in Kustaanhcimo-Sticfcl (KS) coordinates 21, 23] 
originally introduced for the Kepler problem in classi- 
cal celestial mechanics, but also adapted to the hydrogen 



atom [23], transforms the Coulomb potential to a har- 
monic potential with a constraint. In the regularized hy- 
drogen atom the application of the GWP method should 
therefore be capable of yielding exact results when the 
constraint can be handled. The regularization implies 
a fictitious time variable which has been shown to be 
the eccentric anomaly of the corresponding classical orbit 
[2^ . Various approaches have been made to construct co- 
herent states for the hydrogen atom in the fictitious time 
[H [H [13, [H, m, [33, Hlj in analogy with the coher- 
ent states of the harmonic oscillator. These approaches 
construct the coherent states as the eigenstates of the 
lowering operators associated with the harmonic poten- 
tial. 

In this paper we consider the field-free hydrogen atom 
and show that contrary to the dynamics in the real phys- 
ical time the exact propagation of arbitrary initial Gaus- 
sian wave packets in the fictitious time can be described 
analytically and exhibit a nondispersive periodic time de- 
pendence. In SeclUthe Coulomb problem is transformed 
to the problem of the four-dimensional (4D) harmonic os- 
cillator in Kustaanheimo-Stiefel coordinates subject to a 
constraint. The consequences of the constraint for Gaus- 
sian wave packets are discussed in Sec. lIIII and the phys- 
ically allowed set of "restricted Gaussian wave packets" 
is constructed. In Sec. IIVI the analytical time evolution 
is derived for initially three-dimensional (3D) Gaussian 
wave packets in the physical space, and also for two- 
dimensional (2D) and one-dimensional (ID) wave packets 
with cylindrical and spherical symmetry, respectively. 

In the subsequent paper the investigations are ex- 
tended from the pure Coulomb problem to the hydro- 
gen atom in static external electric and magnetic fields. 
Wave packets are propagated by application of the time- 
dependent variational principle in such a way that the 
dynamics is exact for the Coulomb problem and approx- 
imations in the variational approach are only induced by 
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the external fields. Quantum spectra of the nonintegrable 
systems are then obtained by the frequency analysis of 
the time autocorrelation function of the propagated wave 
function. 



Eq. ([7]) is not a standard linear eigenvalue problem. Scal- 
ing the coordinates 



(8) 



II. REGULARIZATION OF THE HYDROGEN 
ATOM 

To make our presentation self-contained we briefly re- 
view the Kustaanhcimo-Stiefel transformation for the 
Coulomb problem. The time-independent Schrodinger 
equation for the hydrogen atom reads 
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(1) 



with A3 the Cartesian form of the Laplacian. A regu- 
larization of the singular Coulomb potential is obtained 
by using Kustaanheimo-Stiefel (KS) coordinates u = 
U2, U3, U4) (2l|. [2^ which are introduced here, ac- 
cording to Ref. [25I, differing by a factor of two from 
the original definition, 



X = U1U3 — U2U4, 

y = uiUi + U2U3, 

1 / 2 I 2 2 2\ 
Z = - Ui + U9 — — 



(2) 



By adding a fourth component with the constant value 
zero to the physical position vector, i.e. x — {x,y,z,0), 
the transformation can be written in matrix notation, 



X = L{u)u , 



with 



L(u) 



U3 -U4 Ui -U2 

1 [ M4 U3 U2 Ml 

2 \ Ui U2 -U3 -M4 

U2 —Ui — U4 U3 



(4) 



The introduction of the auxiliary degree of freedom, 
which renders the originally three-dimensional problem 
four-dimensional, entails a constraint on physically al- 
lowed wave functions ^j, i.e. 



d d d d . , ^ 

U2-r. Ul- M4-- h Uz- V = . 

OUi OU2 OU3 OU4 



(5) 

With r = {x^ +y'^ + z^y^^ = u^/2 the Schrodinger equa- 
tion II]) transformed in Kustaanheimo-Stiefel coordinates 
reads 



■iA4-A,v. = i,^ 



(6) 



where A4 denotes the 4D Cartesian form of the Lapla- 
cian. Multiplication with and reordering of the terms 
yields 



Hi; = ( -iA4 - Eu^ ]iP^2iP 



(7) 



and setting 



E 



1 



(9) 



2n2 ' 

leads to the time-independent Schrodinger equation 

if^= (-^A4 + iu2^ ^ = 271^-, (10) 

which represents the Schrodinger equation of the 4D 
harmonic oscillator subject to the constraint ([5]). The 
scaling parameter, which takes only integer values n = 
1,2,3,..., turns out to be the principal quantum number 
of the hydrogen atom. The 4D isotropic harmonic oscilla- 
tor is invariant under the unitary group U{A) and thus the 
eigenstates of Eq. (jlOp are not unique. The simple prod- 
uct of four eigenstates of the ID harmonic oscillators, i.e., 
the separation of Eq. (flOl) in the Cartesian coordinates 
M2, M3, U4), in general violates the constraint ([5]) and 
therefore represents unphysical solutions. The constraint 
([5]) rather suggests the introduction of two sets of polar 
coordinates, viz. the semiparabolic coordinates 



Ul 

U3 



fj. cos ipf_, 
V cos If I, 



U2 ^ fJ. sm if^ 
— v sin Lpi, , 



(3) with the associated angular momenta {pj = ig^) 



W1P2 - U2V\ 



Eq. (fT^ yields the constraint ([5]) in the form 



(11) 



(12) 



(13) 



The relation between the physical Cartesian coordinates 
and the semiparabolic coordinates is obtained using the 
definitions ^ and (fTT|) 



X — /ii^(cos 93^ cos — sin sin (/3i,) 

— ^iycos{ipfj_ + ifiy) = ^vcosip , 

y — {cos ip sm Lp I, + sin Lp^ cos Lp^) 

= ^vsvaip^ + p^) ~ fivsinp , 



(14) 



with the physical azimuthal angle p — ip^ -\- ip^. In 
the semiparabolic coordinates the Schrodinger equation 
reads 



ip — 2n ijj 



(15) 



3 



with 



1 d d 



1 



p dp dp p2 dtp'^ 



p = p^v 



(16) 



Eq. ()15|) is separated in two uncoupled 2D harmonic os- 
cillators in the coordinates /i, t/s^t and v^Lp^,^ respectively. 
The solution can be taken in the product form 



ij{p,v, ip) = $iV^™(^)$JV„m(j^)e'""^ 



(17) 



where 



1 d d 



1 



'2pdp^dp 2p2 + 2^ 



with p = p,v and Np — 0, 1, 2, 
resentation of the eigenstates is 



(27Vp + |m| + l)$jv,™(p) , (18) 
The coordinate rep- 



m 



7r(7Vp + |m|)' 



pl™l<l(p2)e-i''^ , (19) 



with the associated Laguerre polynomials L^' . For the 
principal quantum number n introduced above we obtain 
the relation 



Nf, + N^ + \m\ + 1 = 1,2,3,. 



(20) 



and therefore via Eq. ^ the correct Rydberg spectrum. 

To perform time-dependent computations it is nec- 
essary to formulate the time-dependent version of the 
Schrodinger equation pO)l . By analogy with the usual 
identification i? — > where t is the physical time, the 
"fictitious time" variable r is introduced, (see, e.g., Refs. 
[l^H^I) as the conjugate variable to the principal quan- 
tum number rt, 



OT 



(21) 



The regularized Schrodinger equation for the hydrogen 
atom in the fictitious time then reads 



i—^ = Hib = ( — A4 + -u^ 
dr ^ I 2 2 



(22) 



In the following the "fictitious time" r will simply be 
denoted as "time" for brevity, whereas t will be named 
"physical time" . 

The form of Eq. ((22| suggests that it can simply be 
solved by Gaussian wave packets in Kustaanhcimo-Stiefel 
coordinates u, i.e. 



5(y, u) 



— g»[(u-q)'4(u-q)+7r-(u-q)+7] 



(23) 



where A is a complex symmetric 4x4 matrix with pos- 
itive definite imaginary part and the momentum tt and 
the center q are real, 4D vectors in the Kustaanheimo- 
Stiefel coordinates. Those vectors represent the expec- 
tation values of the position and the momentum op- 
erator, respectively, i.e. q = {g\\i\g) / {g\g) and tt — 



(.9li^4|5')/(ffl5)- The phase and normalization is given 
by the complex scalar 7. Collectively the parameters 
are denoted by y = (^, 7r,q, 7), which is a set of 
4 X (4 -I- l)/2 -I- 4 -I- 1 = 15 complex parameters when 
the two real vectors q and tt are counted as a single com- 
plex vector. Inserting ()23|) into the Schrodinger equation 
([22|) yields a set of ordinary differential equations for the 
parameters y, which can be solved analytically for any 
given initial GWP. As the potential is harmonic the time- 
dependent GWPs are exact solutions of Eq. . 

The problem, however, is that in Kustaanhcimo-Stiefel 
coordinates the constraint ([5|) on physically allowed wave 
functions must be taken into account. The GWPs in 
general violate that constraint. The question is whether 
the constraint can be fulfilled exactly by a Gaussian at 
all, and, if so, whether GWPs fulfilling the constraint 
are reasonable trial functions in the sense that they still 
present a complete basis set. 

The advantage of the formulation of the Hamiltonian 
in semiparabolic coordinates IjlSp is that the constraint 
([5]) is already incorporated. However, in semiparabolic 
coordinates it is nontrivial to find Gaussian type trial 
functions for the exact solution of the time-dependent 
Schrodinger equation because the Laplacian is not of 
Cartesian form, see Eq. and includes centrifugal 

barriers. In the following we will therefore use the Carte- 
sian type KS coordinates rather than semiparabolic co- 
ordinates to investigate the impact of the constraint ([5]) 
on a Gaussian trial function and to discuss the properties 
of the resulting "restricted Gaussian wave packets" . 



III. RESTRICTED GAUSSIAN WAVE PACKETS 

The regularization of the hydrogen atom in Sec. ITTlhas 
transformed the Coulomb potential to a harmonic poten- 
tial in the Schrodinger equation (fTU]) . The goal now is to 
perform exact wave packet propagation in the fictitious 
time r for the hydrogen atom. The constraint can be 
fulfilled by a 4D GWP if the space of admissible config- 
urations of the parameters y is restricted. The structure 
matrix 



J = 



10 

-10 

-1 

1 0. 



(24) 



allows for the compact notation of the constraint ([S]), 



(25) 



Letting the constraint operator X act on the trial func- 
tion (USD yields 

Xg{y, u) = uJ^[2A(u - q) + tt] .g(y, u) ^ . (26) 

Thus the result is a quadratic polynomial in the coordi- 
nates u, multiplied by the GWP itself. The constraint 
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(l26)l has to be satisfied pointwise for all u € K^. For non- 
trivial wave packets the polynomial in (|26p must vanish, 
and an algebraic equation remains, 



2Aq) = 



(27) 



which is only possible if all coefficients of the second or- 
der polynomial in ([27]) are zero. Let us first investigate 
the term linear in u whose coefficients must vanish, i.e. 
TT — 2A(\. Inserting this condition into the wave func- 
tion ([23]) yields g{y, u) = exp{«[(u - q)^(u + q) + 7]} = 
exp{i[uj4u + 7']} with 7' = 7 — q^q, which means that 
without loss of the variational freedom we can set 



q = , and tt = 



(28) 



because any nonzero q and tt vectors only change the 
scalar 7. The bilinear form in (I27p must also be zero, i.e. 
2uJ^Au = 0. This requires the matrix of the bilinear 
form 



J' A = 



I 



\ 



-ai2 
ail 
ai4 

-ai3 



-022 
ai2 
a24 

-0.23 



-023 —024 
^34 044 

-033 —034 



to be skew-symmetric, i.e., the diagonal elements of J^A 
must vanish, 012 — 034 — 0, and from the off diagonal 
elements we obtain an = a22, 033 = ^44, 124 = ~ai3, 
and 023 = 0-14. With the definitions an = a^, 033 = a^, 
0-13 = ttx, an = ay, the matrix A must be of the form 



A 








ax 


ay 







ay 


~ax 


ax 


ay 


a^ 





\ ay 


-ax 





a^ 



(29) 



Eqs. (|28p and imply that the set of 15 complex pa- 
rameters y of the general GWP ((23)) is reduced to only 
five parameters y — {a^, a^,, ax, ay,"f) for a restricted 
GWP satisfying the constraint ([5])! 

The question arises whether the restricted GWPs form 
a complete basis set, such that any physically allowed 
state can be expanded in this basis. It is not evident 
that a superposition of restricted GWPs whose centers 
are all located at the origin and which only differ by their 
complex widths, is flexible enough to represent arbitrary 
quantum states. The usual form of the resolution of the 
identity [33*] for a continuous basis set of normalized, un- 
restricted GWPs of the form dMl) is 



1 



{2n) 



d7TUqMy)){g{y)\=l 



(30) 



where the width of each GWP basis state is kept fixed. 
Obviously, Eq. ([50]) cannot be applied to the restricted 
GWPs, since both parameters tt and q are set to zero 
(see Eq. ([^5))). However, it is sufficient to require the 
restricted GWPs to be complete in the 3D physical space 
only. To verify the completeness in the 3D space we 



transform the restricted GWP in KS coordinates back 
into the original 3D Cartesian coordinates. 



5(y,x)=e^^"-"+^^ (31a) 

= gi[aM("l+"2) + ''''("3+"4)] 

oi[2ax (til ti3-M2U4)+2ay (ill «4 +1*2 1*3)4-7] (31b) 



X e 



^i{af^fj.'^+a^v^+2a^x+2ayy+'f) 
J[{af^+a,^)r+(a^,-a^)z+2a^x+2ayy+f] 



(31d) 
(31e) 



where we have exploited that in semiparabolic coordi- 
nates is^ = r + z, i''^ = r — z. In Eq. piep the set of 
parameters (a^, ai,,ax,ay) has been replaced by an equiv- 
alent set of complex parameters, defined by 



ai. 



Pr — a^ -r u,u , 

P = iPx,Py,Pz) = i'2ax,2ay,a^ 



(32) 



For Pr — and real- valued parameters Px,Py,Pz the re- 
stricted GWP in Cartesian coordinates (|31ep reduces to 
a plane wave e^^'^. (Note that in that case the imagi- 
nary part of the matrix A is not positive definite and the 
wave function cannot be normalized.) Since plane waves 
are known to form a complete basis it is proved that the 
restricted GWPs, forming a superset of plane waves, are 
also complete, or even over-complete. However, they do 
not form a complete basis set of the 4D harmonic oscil- 
lator pH)) . 

The GWPs dni) satisfy the constraint (0, but "for the 
price" of the condition ([^5]) . i.e., they are localized around 
the origin with zero mean velocity. It seems impossible 
that the time propagation of a single restricted GWP 
exhibits any meaningful dynamics with, in particular, a 
classical limit in the sense of the correspondence prin- 
ciple. However, a wave packet localized in the physical 
coordinate and momentum space of the hydrogen atom 
can be constructed as a superposition of the restricted 
GWPs. The expansion and exact time evolution of wave 
functions in the basis pip , and in addition, in modified 
bases sets for certain symmetry subspaces of the hydro- 
gen atom are the subjects of the next section. 



IV. ANALYTICAL WAVE PACKET DYNAMICS 
IN THE HYDROGEN ATOM 

An arbitrary wave packet of the hydrogen atom can 
be propagated analytically in the fictitious time. This is 
achieved by expanding the initial wave packet in terms 
of the restricted GWPs, whose time-dependence is de- 
rived and shown to be given by simple analytical for- 
mulae. We consider three different cases. In Sec. IIV AI 
the time propagation of Gaussian wave packets in the 
3D physical space is discussed. These wave packets are 
not eigenstates of the angular momentum operator. We 
then introduce the time propagation of wave packets with 
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symmetries, viz. in Sec. II V B] cvlindricallv symmetric 2D 
wave packets with well defined magnetic quantum num- 
ber TO, and in Sec. IIV CI spherically symmetric ID wave 
packets with well defined angular momentum quantum 
numbers I and m. 



A. Propagation of 3D Gaussian wave packets 

The aim is an exact time propagation of arbitrary wave 
functions in the hydrogen atom. The initial wave func- 
tion ■0(0) is taken as a superposition of the restricted 
GWPs ([5T|) . These basis states are then propagated ana- 
lytically in time. We start with the derivation of the time 
evolution of the basis states. Inserting the ansatz (I31ap 
in the time-dependent Schrodinger equation (j22p yields 



d 



(33) 



-7 + itr^ - u{A + 2A^)u - -u^ ) 5(y,x) = . 



The equations of motion for the 4x4 width matrix A 
given in Eq. (j29p and the complex phase factor 7 are 
now directly obtained from Eq. (|34p as 



2 ' 

j = itr A . 



(34a) 
(34b) 



It is important to note that although the parameter set 
of the GWP (j31ap is restricted to only five complex pa- 
rameters (as compared to 15 parameters in (j23p ) the dif- 
ferential equations (I34p still describe the exact dynamics 
of the wave packet without any approximation, i.e., with 
the initial condition the wave packet stays centered 
around the origin for all times, and the width matrix A 
keeps the form ([29)1 because the matrix A in (|34a|) has 
the same structure ([2^ as A itself. 

Both equations can be solved analytically. Equa- 
tion p4ap is solved most easily when two auxiliary com- 
plex matrices B and C are introduced with A = ^BC^^, 
and the initial conditions B{0) = 2A{0) and C(0) = 1 



[34|. Then (|34ap is replaced with the two equations 
C = B and B = -C, or equivalently B = -B. The 
matrices B and C have the same structure as the width 
matrix A with the solution for the matrix B 



Kir) = 
6.(r) = 
byir) - 

and for the matrix C 



c^(r) 
Cxir) 



2a , cos T — sin r 
2a^l cos T — sin r 
2a° cos T , 
2a° cos T , 



COST -f 2a,, sinr , 



(35) 



cos T + 2a° sin t , 
2a° sin r , 

2a°sinr. (36) 
Here and in the following the superscript indicates pa- 
rameters of the initial state at time t = 0. The matrix 
A is obtained from the above definition A = jBC~^ and 
the four elements in the form of Eq. (l32|) read 

Prir) = ^|p^os2r+i[(p°)2-(p°)2-l]sin2r| 



Px{t) 



pI 



Py(j) 



Pz{t) = 



Z{r) 



(37) 



Z(r) ' Z{t)' 

where Z{t) abbreviates the expression 

Z{t) - cos^ T -f [(p°)2 - (p0)2] sin^ t+pI sin 2t . (38) 

With the matrix A at hand it is possible to integrate Eq. 
(|34bp to obtain 7. The quantity e~*'>', i.e., the phase and 
normalization of the wave function, then reads, with the 
initial value 7(0) = 

U{t) = e-*^M = l+[(p°)' - (P°)'] (l-cos2r)+p^sin2r. 

(39) 

The analytical time evolution of the wave function is ob- 
tained by inserting the time-dependent parameters in the 
wave function (|31ep which finally yields 



5(r,y",x) 



■ exp 



2(p" ■ X p°r cos 2t) + [{p°f - (p°)^ - l]r sin 2r 



(40) 



This is an important intermediate result. The time evo- 
lution of a restricted GWP piep has been calculated an- 
alytically and takes the compact form (|40p . The param- 
eters in Eqs. (|37p and (|39p are periodic functions of the 
time T with period tt. This results in a 7r-periodicity of 
the wave function ([40]) . In the physical time, wave pack- 



ets disperse in the hydrogen atom (^,[11]. By contrast, 
the wave packets in the fictitious time show an oscillating 
behavior with no long-time dispersion in t. 



Now that the time evolution of the basis states piep 
is known we can expand an arbitrary initial state in 
this basis. The time evolution of that wave function is 
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then analytically given by the superposition of the time- 
dependent restricted GWPs In general, the expan- 
sion of an arbitrary state V'(x) in an over- complete set 
of Gaussian wave functions g{y,x.) is a nontrivial task. 
A procedure for finding the "optimal" expansion for a 
given number N of basis states is to minimize the devi- 
ation A = IIV'(x) - J2k=i9{y'',^W, e.g., by searching 
for stationary points = with respect to the param- 
eters y'' , k — 1, . . . ,N [36] ■ This procedure presents a 
highly nonlinear minimization problem. Many stationary 
points, i.e., local minima may exist, and the difficulty is 
to find the true global minimum. 

Here, we concentrate on the propagation of 3D Gaus- 
sian wave packets which are localized around a point xo 
with width a in coordinate space and around po in mo- 
mentum space, and present a direct approach for the ex- 
pansion of the 3D GWPs given as 

V'(x) = (27ra2)-3/4 exp {"^^^^^ + * Po • (x - xo) 

(41) 

in the restricted GWPs (I^T]) . The Fourier representation 
of V'(x) reads 



2 \ 3/4 



d^pe 



-cr"(p-po) +ip-(x-xo) 



(42) 



Using the approximation r « x • xo/|xo|, which is valid 
in the vicinity of xq where V'(x) is localized Eq. (j42|l can 
be written as 



V'(x) 



3/4 



X e 



27r3 

*[PrI-+(p-pr^)-x] 

.2 \ 3/4 



(^3pg-o-^(p-po)^-ipxo 



(43) 



2 \ ■5/'' /• 

where the g{y,x.) are the restricted GWPs (|31e[) for 
the set of parameters y given by [pr = const., p — 
PrXo/|xo|,7 — 0). As can be easily shown the restricted 
GWPs (|31ep for constant Pr and 7 = are a complete 
continuous basis as functions of the momentum p, i.e. 



1 



(27r)3 



d'p\9iy)){9iy)\ 



(44) 



Note that the completeness relation ((30|) is valid for 
frozen Gaussians with arbitrary localization in coordinate 
and momentum space while in Eq. (144^ the restricted 
GWPs in KS coordinates are all located around the ori- 
gin but the width matrix A is varied. Note also that Eq. 
(j44|) is exact for pr = 0, i.e., for the restricted GWPs 
given as plane waves, and becomes an approximation for 
Pr 7^ 0. This means that Pr should be chosen as zero or 
close to zero in practical applications. 

In numerical computations it is convenient to approxi- 
mate the 3D Gaussian wave packet (|41l) by a finite num- 
ber of restricted GWPs rather than using the integral 



representation p4)l . This is most efficiently achieved 
by evaluating the integral in (|44p with a Monte Carlo 
method using importance sampling of the momenta with 
a normalized Gaussian weight function 



w{p) = ( — 



3/2 



-a (p~po) 



The initial wave packet then reads 



1 



N 

V'(x) = (2w^)-3/4A^ 

k=l 



.g(y^x) e-'P'-^° 



(45) 



(46) 



with y'^ = {pnp'^ — PrXo/|xo|,0), and the p'^, k = 
1, . . . , N distributed randomly according to the normal- 
ized Gaussian weight function (I45p . The wave function 
■0(x) in Eq. is an approximation to the 3D Gaus- 
sian wave packet (j41l) , and the accuracy depends on how 
many restricted GWPs are included. However, it is im- 
portant to note that a wave packet which is strongly lo- 
calized around xq in coordinate space and around po in 
momentum space can be described even with a rather 
low number N (of the order of 50-100) restricted GWPs. 

The time propagation of an initial state in the fic- 
titious time r is now obtained exactly and fully analyti- 
cally in a simple way. In the unperturbed Coulomb prob- 
lem a superposition of time-dependent restricted GWPs 
is uncoupled, i.e., all basis functions propagate indepen- 
dently. Thus, the time propagation of the initial state 
(|46)) simply reads 

1 ^ 

V'(t,x) = (2w2)-3/4_^^(^^yfe^^) g-.p-^-xo ^ (47) 

fe=l 

with the same parameters y'^ as above and the time prop- 
agation of the restricted GWPs given by Eq. (|40|) . 

Eq. (j47| is the final result of this section, and is il- 
lustrated in Fig. [1] for an initially Gaussian wave packet 
(|4ip whose center is moving in the z — plane. The 
wave packet is expanded and propagated analytically in 
a basis of = 10000 restricted GWPs. This is possible 
since the analytical approach allows for a large number 
of basis states. We note, however, that results of simi- 
lar quality can be obtained using A = 50 ~ 100 GWPs 
only. The probability density in the z = plane is plot- 
ted at equidistant times with step size At — n/b. As 
mentioned above Eq. (|44[) is exact for pr — 0, however, 
in our numerical calculations we choose Pr = ie with a 
small e > 0. The damping enables normalization of the 
restricted GWPs, and improves the convergence of the 
Monte Carlo integral. The initial GWP presented in Fig. 
[T]for T = is centered at xg = (8,0,0) with the mean 
momentum po = (1,2,0), and Pr is set to Pr = ie with 
e = 0.01. Classically the electron with these initial condi- 
tions is running on the Kepler ellipse plotted by dots on 
the bottom of each panel in Fig. [TJ For every time step, 
that part of the ellipse that has been passed by the elec- 
tron so far is shown by a black solid line for time resolved 
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|\lf(x)|^x10'' |\|/(x) 1^x10" \\\i{x)f x^o'^ 




FIG. 1: (Color online) Fictions time propagation of the initially Gaussian wave packet (|41[l located at xo — (8,0,0) with the 
momentum po — (1,2,0) plotted in the plane z — 0. The classical Kepler ellipse with the same initial conditions is indicated 
by the dotted curve on the bottom of each panel. For time resolved comparison, that part of the ellipse that has been traversed 
by the particle so far in each plot is shown by a solid black line. Although the Gaussian wave packet does not stay Gaussian 
during the period it follows in general the classical path and is recovered after one period t = tt, indicating the periodicity of 
the wave packet. 



comparison. The position of the maximum of the prob- 
ability density agrees well with the classical position of 
the electron on the ellipse for all times. The 7r-periodicity 
of the motion is reflected by the coincidence of the wave 
packet after one period at r = tt with the initial GWP 
at r = 0. 

B. Propagation of 2D cylindrically symmetric wave 
packets 

In this section basis functions based on the restricted 
GWP (j31cp with a well defined angular momentum com- 
ponent Iz = m are derived. This case is especially impor- 
tant when a cylindrically symmetric external field, e.g. a 
magnetic field, is applied to the hydrogen atom, as dis- 
cussed in the following paper [s^]. First wave packets 
with definite 1^ are constructed and their exact, analyt- 
ical dynamics in the hydrogen atom is discussed. Then 
we introduce a procedure to expand quantum states of 
defined Iz in terms of the basis states. 



A 2D cylindrically symmetric restricted GWP (|3ip is 
obtained by setting a^; = Oj, = in (I31cp and introducing 
parabolic coordinates ^ = ^'^^r + z,ri = iy^=r~z, i.e. 

go{y, x) = e'(°"^'+'''^ = g*(p«C+P.')+7) ^ (43) 

with the parabolic momenta = and p^i — a^. These 
states are axisymmetric, and have the quantum number 
TO = 0, but can be generalized to arbitrary m by setting 

g„(y,x) = (;x:.)l™le'(°"'''+'''^'''+'')e™'^ 

= (^77)l'"l/^e*^''^''^''''''^''-' e'™'^ . (49) 

The time-dependent parameters are y = (a^,a^,7) or 
equivalently y = (j>(^,Pri,j)- The quantum number m is 
constant. As will be shown, the wave packet still 
presents an exact solution of the regularized Schrodinger 
equation. The Laplacian in semiparabolic coordinates 
(jl6p and the time derivative acting on the wave packet 
(pg]) yield 
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d 

i^3m(y,x) = (-d^^^ - ct^i^^ - 7) .g™(y,x) , (50) 

and thus the time-dependent Schrodinger equation reads [(d^ + 2a^^^ + i)^^ + (d^ + 2a^i^^ + + 7 — 2i(a^ + 
a„){l + |m|)] (7m(y,x) = 0. This equation is solved exactly if the Gaussian parameters obey the equations of motion 

= -2a' - \ , (51a) 

d, = -2a2 _ i , (51b) 

7 = 2i(a^ + a^)(l + |m|) , (51c) 
or using the matrix notation (|29p for v4 (with = ay = 0) 

i = -2A^ - il , (52a) 
7 = itrA(l + |m|) . (52b) 

The equations of motion for the two nonzero complex width parameters and au remain completely unchanged as 
compared to the restricted GWP in Sec. IIV Al The only change is the additional factor of (1 + \m\) in Eq. (j52bp for 
the phase parameter 7. The solution of Eq. (|52ap is 

= ^{2(a°-al^) + 2(a;i + a;^)cos2r-(l-4aX)sin2r} , 

aAr) = ^{2(aO-aO) + 2(aO +aO)cos2T-(l-4aOaO)sin2r} , (53) 

with 

Z{t) = 2 [1 + Aa^al + (l - 4a51a°) cos 2r + 2 (a° + a^) sin 2t] . (54) 

The solution of Eq. (j52bp is the solution of Eq. (j34bp multiplied by the factor (l + |m|), and the phase and normalization 
factor of the wave function reads, with 7(0) — 

■ -. s \rn\ + l 

. (55) 
The time evolution of the wave packet P5|) then is given by 

ff„(T, y«, I.) = L__(^j.)l™l e»[a.(r)M=+a.(r).^] ^^m^ (56) 

The time-dependent basis states (1551) are the analogue of the restricted GWPs (UOl) for the propagation of wave packets 
with constant magnetic quantum number m. They obey the constraint ([5]), but they are not sufhciently general to 
describe the dynamics of 2D wave packets localized around a given point (^qj VojPe.oyPvo) parabolic coordinate 

phase space. Such localized states can now be constructed in a similar way as described in Sec. IIV Al We use the 
formal plane wave expansion of a Gaussian wave packet in parabolic coordinates, viz. 

ipQ{^,r]) = Ac^j>< ^«PCo(C- Co) +«P7,o(»7- 

A(T 



TT 

2 



dp^ J dpy^er'^ liP(-PSa) +iP'i-Pva) ]-iiPiio+PvVo) giiPii+PriV) 
— J dp^ J dprj g-o'^[(P«+"-P«o)^ + (P'J+*^~P'ra)^l+^(?o+''o)~*(Pe'5o+P'j''o) gi[(Pi+i<^)^+(Pr,+i<^)v)] 



Ao'' 
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where ^ is a normalization factor, e has been introduced (without any approximation) as an additional free parameter, 
and go{y, x) is the cylindrically symmetric restricted GWP for the set of parameters y = {p^ + «e,Pr; + 7 = 0). 
A value of e > guarantees that go{y,x) can be normalized. From ([57| an initial state with given magnetic quantum 
number m is obtained as 

V'™(x)-(e7?)l'"l/2^o(C,'?)e™^. (58) 

In numerical computations the integrals in (|57p are approximated employing a Monte Carlo technique in the same 
way as explained in Sec. IIV Al We obtain 

.2 -1 ^ 

TT N 

k=l 

.21 N 

= ^ y^ffm(y'°,x) c~''''^"^°+P^>''"^+"'^°+"°^~^''^"''^''^~'°''c^ + 'P^~'°"'3^^+^'^''"" , (59) 

with sampling points p^, pj^ randomly distributed around p^g, Prjo according to the weight function w{p) = 

{a / ^/n)e^°' ^p^po) . Finally, the replacement of the initial basis states g,„(y'^,x) with the corresponding time- 
dependent solutions ([55)1 yields 

^fT,x.) = g„(r,y^x)e-'(P?^«+P^«)+^(«°+"°)-2^"'^[(pt-P«o)+K-P'-o)]+2-'^' , (60) 

TT iV 

fc=l 



with the parameter sets y'^' = (pt + ie,p'i + ie, 7 = 0). 



In Fig. [2] the expansion of a GWP (I57|) is shown with 
the center at Co = Vo = 25.0 and the momenta p^^ — 
0.535 and p^jg = —0.117 or in terms of cylindrical coor- 
dinates po = 25.0, zq = 0.0 and ppg = 0.419, = 0.652 
and the width a = 4.472. Note that a Gaussian shape of 
a wave packet in parabolic coordinates is nearly Gaus- 
sian also in cylindrical coordinates (see e.g. the wave 
packet at r = in Fig. [21). The wave function shown 
has zero angular momentum component 1^ = 0. For rea- 
sons of presentation the originally positive radial coordi- 
nate p is extended to negative values, and the symmetry 
ip{—p) = ip{p) is used. The probability density p\tp{p, z)P 
is plotted in the (p, z) plane. The Kepler ellipses plotted 
on the bottom in each panel of Fig. [5] show the corre- 
sponding classical motion of the particle with the initial 
conditions given above. The second ellipse again is ob- 
tained by reflection symmetry as the intersection of the 
torus, which is obtained from rotating the ellipse around 
the z axis. At times r w 27r/5 and t w 47r/5 the high 
probability density close to the z axis leads to interfer- 
ence patterns. After one period t — n the initial wave 
function r = is recovered. A number of = 5000 
modified basis states (|49l) with e ~ 0.05 are employed. 
Results for the case m 7^ are not shown since they dif- 
fer only qualitatively by avoiding the crossing of the z 
axis due to the rotational barrier. 



C. Propagation of ID spherically symmetric wave 
paclcets 

The procedure of the two previous subsections is ap- 
plied to quantum states with conserved angular momen- 
tum. First an extension of the basis states piep to ba- 
sis states with well defined angular momentum quantum 
numbers Im is presented, and they are shown to be exact 
solutions of the time-dependent Schrodinger equation of 
the regularized hydrogen atom. Then the procedure of 
expanding states with definite Im in the constructed ba- 
sis states together with an example are presented. For 
radial symmetry the complex width matrix A (1^^)) of the 
restricted GWP must be a multiple of the identity ma- 
trix A = al, i.e. = CLy — 0, a^i = = a = Pr/'i- The 
restricted GWP reduces to 

gm{r) = e'(2"''+'') = e'^^-'^+T) . (61) 

This is a suitable basis state with vanishing angular mo- 
mentum. The correct extension to arbitrary angular mo- 
menta is given by 

gim{r, e, ip) = r'e*(2a.+7)y^^(0^ ^) ^ (62) 

where Yim{0,y^) denotes the spherical harmonics. In- 
sertion of the ansatz (|62)) into the time-dependent 
Schrodinger equation pop in spherical coordinates 

d / 92 p \ 

i-g^giTnir,d,ip) = ( -^»' + ~ + ^) 9imir,0,ip) (63) 

yields 

-t- l)a + 7 -I- r(l -I- (2a)2 -|- 2a)] .g,„,(r, 9,ip) = 0. 

(64) 
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FIG. 2: (Color online) Fictitious time propagation of a GWP (|57p with magnetic quantum number m — 0. The part on the 
negative p axis is obtained by reflecting the positive part at the z axis. The wave packet runs along the classical Kepler ellipse 
with the corresponding initial values. For details see text. 

The basis sates present an exact solution of the 
Schrodinger equation provided the time-dependent pa- 
rameters obey the equations of motion 

a = -2a2-i, (65a) 
j^Aia{l + l), (65b) 

with the analytic solution 

4a0 cos 2t - (1 - 4(a°)2) sin 2t 
~ 2 [1 + 4(a0)2 + (1 - 4(a0)2) cos 2t + sin 2t] 

(66) 

and 

e-'''(^) (67) 
= [1 + 4(a")^ + (1 - 4(a°)2) cos 2t + 4a° sin 2t] /2 . 

Wave packets with well defined angular momentum quan- 
tum numbers I and m can now be expanded in the basis 
(15^ . For the radial coordinate r the same procedure as 
introduced in Sec. IIVBI for the parabolic coordinates ^ 
and T] can be applied. The plane wave expansion of a 



Gaussian wave packet localized around (ro,Prn) reads 
V-ooW = ^expj- ^ +ipr„{r-ro)y (68) 

= / dp.e-'^'(^"-+-^"-o)'+-«-^'"-°goo(y,r) , 

where 500 (Yi^) is the spherically symmetric restricted 
GWP for the set of parameters y = (pr + 7 = 0). 
An initial state with given angular momentum quantum 
numbers I and m is obtained as 

=r'^ooir)Yi„,{9,^) . (69) 

Using the Monte Carlo evaluation of the integral in ([69]) 
and the replacement of the initial basis states (|6ip with 
the corresponding time-dependent solutions 500 (t, = 
exp{i[2a(T)r -|- 7(r)]} (with the time-dependent parame- 
ters given in Eqs. (|66p and ((68)) ) we finally obtain 

Act 1 

i^im{T,x) = ^5/,„(r,y'',x) 

^''^ k=l 

X g-jPr>o+ei-o-2io-2e(p^-pr(,)+(T^e^ ^ /yQ\ 
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FIG. 3: (Color online) Fictitious time propagation of the 
wave function (|69|l with ip{r) given by the Gaussian (|69|l with 
ro = 10.0, Pro ~ —0.5 and a = 3.0 for angular momentum 
quantum numbers (a) Z = 0, m = 0, and (b) I = 5, m — 0. 
An interference pattern is observed at the turning points at 
approximately r = 2n/5 in both panels. The initial wave 
function is expanded in A'^ = 10000 basis states (|62|l with 
e = 0.2. 

with sampling points randomly distributed around 
Pro according to the weight function w{p) — 
{a/y/n)e~'^^'-P^P°^' , and the parameter sets y'' — (pj? + 
ie,7 = 0). 

The results of the propagation of the wave function 
ipim{r, 0, ip) = r'^ i>oo ir)Yim{0, (fi) with ^poo{r) given in Eq. 
([69)1 and the initial values tq = 10, Pro = —0.5 and the 
width (T = 3 is presented in Fig. [3] for different times 
< r < TT. The imaginary parts of pj5 = 2a'' are set to 
e ~ 0.2 and the number of basis states ([5^ is iV = 10000. 
In Fig. [3l^a) the angular momentum is set to to = Z = 
and in Fig. EJ^b) the components of the angular momen- 



tum are / = 5 and to = 0. Due to the negative initial 
value of the radial momentum p^g the wave is initially 
running towards the nucleus located at the origin. The 
wave function with zero angular momentum in Fig. [2I^a) 
comes close to the origin r — 0. Similar to the radial 
symmetric case in Sec. IIV SI there appears to occur some 
interference pattern due to the overlapping parts of the 
incoming wave function at the inner turning point (see 
T = 27r/5). In the nonvanishing angular momentum case 
(Fig. ^h)) the barrier of rotational energy prevents the 
wave function from reaching the nucleus. Instead there is 
a turning point whose distance from the nucleus increases 
with growing angular momentum. Again at r = 27r/5 an 
interference pattern is observed close to the inner turning 
point. In both panels the maximum of the probability 
density overshoots the position of the initial maximum 
To = 10 at T = 47r/5 due to the initial kinetic energy and 
returns to the initial wave packet after the period t — tt, 
indicating the periodicity of the wave function. 



V. CONCLUSION 

In this paper we have derived the wave packet dy- 
namics for the field-free hydrogen atom in a fictitious 
time variable. The Coulomb problem has been trans- 
formed to the four-dimensional harmonic oscillator in 
Kustaanheimo-Stiefel coordinates with a constraint. The 
"restricted Gaussian wave packets" obeying that condi- 
tion have been constructed and their exact time depen- 
dence is calculated analytically. The wave packets with 
and without symmetries exhibit a nondispersive periodic 
behavior in the fictitious time. 

It should be noted that the wave packet propagation in 
the fictitious time substantially differs from the physical 
time dynamics, and thus cannot provide the analytical 
propagation of dispersive wave packets in the physical 
time [3, [l^. Nevertheless, the fictitious time dynam- 
ics can be used to solve the Schrodinger equation for 
Coulomb systems with strong time-independent pertur- 
bations, e.g., the hydrogen atom in static external electric 
and magnetic fields. The restricted Gaussian wave pack- 
ets are the basis for the application of the time-dependent 
variational principle to the hydrogen atom in external 
fields and the computation of quantum spectra by fre- 
quency analysis of the time autocorrelation function in 
the following paper [s^l • As a consequence of using the 
fictitious time variable the method is exact for the field- 
free hydrogen atom and approximations in the variational 
approach are only induced by the external fields. 
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